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Sir: 

Prior to action on the merits, using marked-up versions of the Specification and 
Claims and also clean copies of such, we respectfully request consideration of the 
following amendments and remarks: 

In the Specification (Marked-Up Version) 

On page 1, before the first paragraph, kindly insert the following heading: 
FIELD OF THE INVENTION 

On page 1, between the first and second paragraphs, kindly insert the following 
heading: 

BACKGROUND 

Please replace the first full paragraph on page 2 with the following: 
E. W. Cambel et al. in "J. Physiol.", Vol. 197, 1968, p. 551 and R. De Valois et 
al. in "Vision Res.", Vol. 22, 1982, p. 545, describe psychophysiological experiments 
demonstrating that the human visual system breaks the retinal image down into a number 
of filtered images, each of which [containing] contains intensity variations spanning a 



single local frequency or orientation range. Feature extraction procedures were developed 
on this basis using a multi-channel filtering with so-called Gabor filters (see A.K. Jain et 
al. in "Pattern Recognition", Vol. 24, 1991, p. 1167 ff). In these procedures, a window 
with a specific filter function is incrementally pushed over the image, and the filter 
response is further evaluated as the local average for the respectively examined window 
area. The disadvantage to these procedures is that a linear feature filtering takes place, 
which enables only limited pattern recognition. For example, a linear filter does not make 
it possible to recognize texture differences within an image according to Fig. 7a (see 
below). Filter procedures in which respectively determined averages for a window area 
are subjected to nonlinear further processing were developed to counter this problem. 
However, this requires additional information about the patterns to be recognized to 
achieve a sufficient reliability during texture recognition. 

Please replace the first paragraph on page 3 with the following: 
In addition to acquiring the image features according to one of the available 
procedures, there is another problem involving the evaluation of information contained in 
the image features to achieve the desired image segmentation. To this end, a so-called 
"cluster analysis" is usually performed in a feature space, in which each axis represents 
one of the examined features. A characteristic marker or label is allocated to each image 
point, wherein image points with the same label belong to the same feature class or 
structure, while image points with varying labels must be allocated different structure. 
The available cluster analyses are described, for example, by R.C. Dubes in "Handbook 
of pattern recognition and computer vision", published by C.H. Cheng et al., World 
Scientific Publishing 1993, p. 3 ff., and B.D. Ripley in "Pattern Recognition and Neural 
Networks", Cambridge University Press, 1996. Unmonitored cluster algorithms, i.e., 
cluster algorithms without supervision, initially process unmarked data, and therefore 



require that the following two problems be solved. 

Please replace the first full paragraph on page 4 with the following: 
DE-OS 43 17 746 discloses a spatial filter procedure for recognizing structures in 
n-dimensional images based on the concept of the so-called " isotropic scaling factor a." 
Scaling factor a describes the change in point density (gradient) around an examined image 
point by indicating the surrounding point number as a function of the distance from the 
examined image point. DE-PS 196 33 63 describes the expansion of this spatial filter 
procedure to recognizing the orientation of structures in n-dimensional images. This 
concept introduces anisotropic scaling factors a jS that are characteristic for the point density 
variation after projection on specific directions in space. Spatial filtering based on the 
scaling factors constitutes a nonlinear procedure, which has been successfully used during 
pattern recognition and also nonlinear time sequence analysis. However, this did not 
enable image segmentation relative to different textures in the examined area. 
Please replace the paragraph spanning pages 4 and 5 with the following: 
[The object of the invention is] Accordingly, it would be advantageous to provide 
an improved method for segmenting a point distribution relative to textures that overcomes 
the mentioned disadvantages to conventional procedures, and in particular has a high 
sensitivity and reliability, requires as little advance information as possible about the point 
distribution, and can be used on as wide a scale as possible in the most varied of tasks, 
both in terms of conventional optical images, and of lower or higher dimensional 
formations. [The object of the invention is] It would also be advantageous to provide a 
device for realizing such a method and applications of the method. 
Please delete the first full paragraph on page 5. 

On page 5, before the second full paragraph, kindly insert the following heading: 
SUMMARY OF THE INVENTION 



Please replace the second paragraph on page 8 with the following: 

In general, the following properties can be allocated to the textures. A local 
sequence repeats over an image region, which is large compared to the size of the local 
sequence. In this case, the sequence consists of an ordered arrangement of basic 
constituents, which form roughly equal units as so-called ^micro-patterns, " have the same 
size within the texture, and can be characterized by specific local properties. 

On page 9, before the first full paragraph, please insert the following heading: 
BRIEF DESCRIPTION OF THE DRAWINGS 

Before the paragraph spanning pages 10 and 11, please insert the following 
heading: 

DETAILED DESCRIPTION 



In the Claims (Marked-Up Version) 
1. (Amended) A method for segmenting a point distribution with numerous 
points into partial areas, which each exhibit specific structural elements[, with the 
following steps] comprising : 

— » — ► 

a) determining for each point (ft), a feature vector (x) [is determined,] 
whose components are determined based on at least several scaling factors 
belonging to [the] a respective point; 

b) [the] determining accompanying feature vectors (p t ) [are determined] 

for a predetermined number of reference points (x* ) of the point distribution 
for which [the] allocation to one of the structural elements is given, and 
texture classes each corresponding to the underlying structural elements are 
formed out of the feature vectors of the reference points; 

c) [for all remaining points (p") of the point distribution that are not 
reference points,] determining the distance between the respective point and 

each of the texture classes [is determined] for all remaining points (p") of the 
point distribution that are not reference points based on a distance measure 
in the feature space, which is fixed by the components of the feature vectors; 

d) allocating each of the points (p") [are each allocated] to [the] a texture 
class for which [the] a lowest distance was determined; and 

e) [the] forming partial areas of segmentation [are formed] out of [the] 
respective reference points belonging to a texture class and the points 
allocated in step d). 



2. (Amended) [A] The method according to claim 1 , in which several isotropic 
and anisotropic scaling factors (a) are determined as components of the feature vectors 

ft). 

3. (Amended) [A] The method according to claim 1 , in which expected values 

of several scaling factors are determined as components of the feature vectors (Xj). 

4. (Amended) [A] The method according to claim 2 [or 3], [in which the] 
further comprising determining anisotropic scaling factors ascertained for a point [are 
determined] in various coordinate systems rotated relative to each other. 

5 . (Amended) [A method according to claim 1 , in which an ellipsoidal distance 
measure is used as a local distance measure in step c) for each texture class] The method 
according to claim 3. further comprising determining anisotropic scaling factors 
ascertained for a point in various coordinate systems rotated relative to each other . 

6. (Amended) [A method according to claim 1, in which a Euclidian distance 
measure is used as a shared distance measure in step c) for all texture classes] The method 
according to claim 1 , in which an ellipsoidal distance measure is used as a local distance 
measure in step c) for each texture class . 

7. (Amended) [A method according to claim 1 , in which the segmented partial 
areas are displayed, temporarily stored and/or further processed] The method according 
to claim L in which a Euclidian distance measure is used as a shared distance measure in 
step c) for all texture classes . 

8. (Amended) [A method according to claim 7, in which the size of the partial 
areas is quantitatively acquired simultaneous to displaying the segmented partial areas] The 
method according to claim 1. further comprising displaying, temporarily storing and/or 
further processing segmented partial areas . 



9. (Amended) [A image segmentation arrangement for segmenting a point 
distribution out numerous points, with a measuring device (1) for acquiring the point 
distribution, a filtering device (2) for scanning and filtering the point distribution, an input 
device (3), a calculating device (4) and an output device (5), wherein the filtering, input 
and calculating devices (2, 3, 4) are designed to execute a method according to one of 
claims 1 to 8] The method according to claim 8. further comprising quantitatively 
acquiring the size of the partial areas substantially simultaneously to displaying the 
segmented partial areas . 

10. (Amended) [Method] The method of [using a device] according to claim [9] 
1. which processes [to process:] 

images of medical or biological objects; 
images of materials; 

point distributions of complex static systems; 

point distributions that reproduce system status of complex, dynamic systems ; 

and 

time patterns for dynamic systems. 

11. (New) An image segmentation arrangement for segmenting a point 
distribution out of numerous points, with a measuring device for acquiring the point 
distribution, a filtering device for scanning and filtering the point distribution, an input 
device, a calculating device and an output device, wherein the filtering, input and 
calculating devices execute a method according to claim 1 . 



Remarks 

We have amended the Specification and the claims to place them into better form 
for examination on the merits and allowance. 

Passage to the appropriate art unit for examination on the merits is respectfully 
requested. 
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METHOD AND DEVICE FOR SEGMENTING A POINT DISTRIBUTION 



The invention relates to a method for segmenting a point 
distribution into partial areas with varying structural 
properties, and a device for executing such a segmentation 
procedure . 

Image segmentation, i.e., dividing an image into segments 
or partial areas on the basis of specific image features 
shared within a partial area, is one of the most important 
tasks in image processing technology. In the simplest case, 
image segmentation is based on the recognition of 
differences in gray scale values, e.g., within the 
environment of a considered image point, or on edge 
detection techniques. However, only simply structured 
images with two-dimensionally extended, homogeneous image 
elements can be segmented in this way. However, practical 
tasks, e.g., image processing in medical or materials 
technology, involve more complex image structures, e.g., in 
the form of differentiated, repeating gray scale value 
patterns or indistinct delineations of image elements, 
which cannot be acquired using the simple segmentation 
techniques . 

A series of feature extraction procedures were developed to 
process more complex images (see M. Tuceryan et al in 
^Handbook of pattern recognition and computer vision'", 
published by C.H. Cheng et al . , World Scientific 
Publishing, 1993, pp. 235 f f ) , in which local features 
corresponding to different image structures or textures are 
first extracted, and a feature vector composed of these 
features is then allocated to each image point. Cluster 
techniques are then used in the feature space to divide the 
feature vectors of the image points into a few classes, 
after which the image is segmented into corresponding 
partial areas on this basis. These procedures are not 
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parameter-f ree, so that a high level of effort must in some 
cases be expended on preparation for specific image 
segmentation in order to determine optimal input 
parameters. Statistical analyses are also known for 
evaluating local gray scale value distributions of the 
image structures (see R.M. Haralick in ^Proceedings of the 
IEEE", Vol. 67, 1979, p. 786 f f ) , wherein correlation 
techniques are used for structure acquisition, for example. 
These conventional techniques are associated with 
disadvantages relating to data processing outlay and 
reliability. 

E.W. Cambel et al. in „J. Physiol.", Vol. 197, 1968, p. 551 
and R. De Valois et al. in „Vision Res.", Vol. 22, 1982, p. 
545, describe psychophysiological experiments demonstrating 
that the human visual system breaks the retinal image down 
into a number of filtered images, each of which containing 
intensity variations spanning a single local frequency or 
orientation range. Feature extraction procedures were 
developed on this basis using a multi-channel filtering 
with so-called Gabor filters (see A.K. Jain et al. in 
„Pattern Recognition", Vol. 24, 1991, p. 1167 ff ) . In these 
procedures, a window with a specific filter function is 
incrementally pushed over the image, and the filter 
response is further evaluated as the local average for the 
respectively examined window area. The disadvantage to 
these procedures is that a linear feature filtering takes 
place, which enables only limited pattern recognition. For 
example, a linear filter does not make it possible to 
recognize texture differences within an image according to 
Fig. 7a (see below) . Filter procedures in which 
respectively determined averages for a window area are 
subjected to nonlinear further processing were developed to 
counter this problem. However, this requires additional 
information about the patterns to be recognized to achieve 
a sufficient reliability during texture recognition. 



In addition to acquiring the image features according to 
one of the available procedures, there is another problem 
involving the evaluation of information contained in the 
image features to achieve the desired image segmentation. 
To this end, a so-called cluster analysis is usually 
performed in a feature space, in which each axis represents 
one of the examined features. A characteristic marker or 
label is allocated to each image point, wherein image 
points with the same label belong to the same feature class 
or structure, while image points with varying labels must 
be allocated different structures. The available cluster 
analyses are described, for example, by R.C. Dubes in 
^Handbook of pattern recognition and computer vision", 
published by C.H. Cheng et al., World Scientific Publishing 
1993, p. 3 ff., and B.D. Ripley in ^Pattern Recognition and 
Neural Networks", Cambridge University Press, 1996. 
Unmonitored cluster algorithms, i.e., cluster algorithms 
without supervision, initially process unmarked data, and 
therefore require that the following two problems be 
solved. 

First, it is important with respect to the reliability of 
image segmentation that the correct cluster number be 
selected and confirmed. Second, it must be ensured that the 
labels provided by the algorithm reflect physically 
meaningful features. Another disadvantage to the 
unmonitored procedures is that they are based on minimizing 
an image-encompassing (global) energy function through 
iterative methods. This leads to a tendency toward 
solutions that yield balanced cluster occupations (see A.M. 
Bensaid et al. in ^Pattern Recognition", Vol. 29, 1996, p. 
859 ff ) . 

The problems discussed above are encountered not just 
during the segmentation of images that represent the 
optical mapping of a scene with material objects. Instead 
of the two-dimensional formation examined in the process, 
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in which a gray scale value is allocated to each image 
point, for example, an image in the broadest sense can also 
be a lower or higher dimensional formation, in which each 
image point is first defined by a number of coordinates 
corresponding to the dimension of the formation, and a 
specific number of image features (measured values) is 
allocated to each image point. In addition to spatial and 
time dimensions, the formation dimensions can be expanded 
to include any additional feature axes. In the broadest 
sense, then, the examined systems encompass all physical, 
y, chemical or biomedical processes or materials whose state 

O or individual features can be characterized with a set of n 

H parameters according to the dimension. The systems can be 

43 invariable (static) or chronologically variable (dynamic) 

j~J within the examination period. In the latter case, the time 

is one of the n parameters. 

r; DE-OS 43 17 746 discloses a spatial filter procedure for 

fy recognizing structures in n-dimensional images based on the 
'*t concept of the so-called isotropic scaling factor a. 

fy Scaling factor a describes the change in point density 

(gradient) around an examined image point by indicating the 
surrounding point number as a function of the distance from 
the examined image point. DE-PS 196 33 63 describes the 
expansion of this spatial filter procedure to recognizing 
the orientation of structures in n-dimensional images. This 
concept introduces anisotropic scaling factors a j± that are 
characteristic for the point density variation after 
projection on specific directions in space. Spatial 
filtering based on the scaling factors constitutes a 
nonlinear procedure, which has been successfully used 
during pattern recognition and also nonlinear time sequence 
analysis. However, this did not enable image segmentation 
relative to different textures in the examined image. 

The object of the invention is to provide an improved 
method for segmenting a point distribution relative to 



textures that overcomes the mentioned disadvantages to 
conventional procedures, and in particular has a high 
sensitivity and reliability, requires as little advance 
information as possible about the point distribution, and 
can be used on as wide a scale as possible in the most 
varied of tasks, both in terms of conventional optical 
images, and of lower or higher dimensional formations. The 
object of the invention is also to provide a device for 
realizing such a method and applications of the method. 

This object is achieved by means of a method or device with 
the features set forth in claims 1 and 9. Advantageous 
embodiments and applications of the invention are described 
in the subclaims. 

in a first important aspect of the invention, a method for 
the partially monitored segmentation of point distributions 
is provided in which structure or texture features are 
determined for each point using the mentioned concepts of 
isotropic and anisotropic scaling factors, and a cluster 
procedure with partial supervision is used for 
segmentation. In the cluster procedure, each point is 
allocated to one of the predetermined texture classes, 
proceeding from a predetermined number of known classes of 
structural elements (texture classes) and their allocation 
to specific points (reference points, points with labels) 
and a distance measure that defines the difference between 
the respective texture features and each of the texture 
classes for each of the remaining points (points without 
labels) . The number of reference points is selected as a 
function of application, in particular as a function of 
image size, in such a way that enough points can be taken 
into account for a statistically reliable assessment for 
each texture class. Each texture class should preferably 
encompass approx. 30 to 40 (or even 100) reference points. 



- 6 - 

The mentioned allocation takes place for the initially 
unclassified points without labels via the evaluation of 
the distance measure, preferably by allocating the 
respective point to the texture class from which it 
exhibits the lowest distance. The cluster procedure 
realized according to the invention is referred to as 
partially monitored or as a procedure with partial 
supervision, since the limited number of reference points 
are classified with a known texture allocation. This 
ensures that image segmentation proceeds from a sufficient 
cluster number and physically meaningful labels. One 
special advantage to the invention lies in the fact that 
image segmentation is highly reliable, even if the number 
of reference points is substantially lower (the percentage 
of reference points can lie under 1 %, e.g., 0.1 %) than 
the number of label-free points. 

In a preferred embodiment of the invention, the distance 
measure is defined for each texture class specifically as a 
function of the orientation and form of the point set that 
was allocated to the predetermined texture within the 
framework of partial supervision. However, use can also be 
made of more easily defined distance measures defined 
jointly for all texture classes in the global feature range 
of the point distribution. 

All points in the point distribution that are allocated to 
a texture class form a texture segment, which is 
subsequently displayed or subjected to further processing. 

In another important aspect of the invention, a device is 
described for realizing the texture segmentation procedure. 
This arrangement encompasses a device for measuring the 
point distribution and the feature of the respective system 
state belonging to each point, a filter device with means 
for scanning the examined points in the point distribution, 
means for counting points in the environment of examined 



points, means for acquiring predetermined scaling factors, 
and means for statistically processing the scaling factors, 
an input device designed to allocate predetermined 
reference points for which a texture class affiliation is 
known to the corresponding texture classes, a calculating 
device for ascertaining and evaluating distance measures 
for the texture features of the remaining points relative 
to the texture classes, and an output device, with which 
the texture segments are displayed, temporarily stored or 
passed on for further processing. 

In the broadest sense, the point distributions processed 
according to the invention can represent system states in 
an n-dimensional state space. The point distribution 
represents a two- or higher dimensional mapping of the 
system state, so that reference is generally made to image 
segmentation and image points below. The segmented images 
can also encompass signal or amplitude progressions as a 
function of a reference parameter (e.g., time, energy, 
etc.) or optical gray scale value and/or color images. In 
addition to image patterns, the examined systems can also 
include in particular materials, mechanical devices or 
biological systems. Depending on the application, a system 
state is acquired via measures involving actor and sensor 
technology, analysis and recording or signaling. The 
potentially required actor technology involves measures for 
generating system reactions representative for 
characteristic states, e.g., the excitation of mechanical 
oscillations in an object under examination, or the 
initiation of evoked potentials in neurological systems. 
The sensor technology encompasses the detection of system 
features with respect to the relevant n parameters, and the 
display of features in a high dimensional feature space, 
e.g., via the storage of suitable value groups that are 
allocated to the features. 



In the examined point distribution of image, a complex, yet 
definable image structure is called a texture. A texture 
forms an image area or image region, in which the image 
structure can be traced back to a recurring pattern that 
accommodates elements according to a configuration rule. In 
other words, a specific (constant) texture can be allocated 
to an image region if a set of local statistics or other 
local properties of the image features is constant or only 
changes a little. An image texture is described by the 
number and type of its (gray) tone basic elements, and the 
spatial arrangement of these elements. 

In general, the following properties can be allocated to 
the textures. A local sequence repeats over an image 
region, which is large compared to the size of the local 
sequence. In this case, the sequence consists of an ordered 
arrangement of basic constituents, which form roughly equal 
units as so-called micro-patterns, have the same size 
within the texture, and can be characterized by specific 
local properties. 

In applications involving image processing, e.g., in 
medical or material technology, the recognition and 
segmentation of textures makes it possible to differentiate 
between various object areas. For example, when examining 
ultrasound images of the inner tissues, texture 
segmentation makes it possible to differentiate between 
tumors and healthy tissue, and determine the size of the 
tumor tissue. In this case, the automation of the 
recognition and size indication process is of particular 
interest. Texture segmentation according to the invention 
generally allows both the allocation of textures to 
specific points and a determination of the size of the 
respective image regions with a constant texture. 

The invention offers the following advantages. The texture 
segmentation procedure is parameter-free and non-iterative, 
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as explained in detail below. No free parameters need to be 
optimized during the classification of image features. As a 
result, the method has a high segmentation rate and 
reliability. The segmentation result is slightly sensitive 
to the specific selection of boundary conditions during 
image evaluation. A nonlinear filter technique for texture 
feature extraction is introduced for the first time. The 
segmentation procedure can be used without any problem for 
any tasks by adjusting the nonlinear filtering process as a 
function of application. For the first time, textures can 
be simultaneously recognized and quantitatively analyzed in 
relation to their expansion. 

Other advantages and details of the invention will be 
described below drawing reference to the attached 
flowcharts and drawings. Shown on: 

Fig. 1 is an example for a pixel image with four natural 
textures (a) and the result of texture segmentation 
into four feature classes (b) , 

Fig. 2 is a flowchart to illustrate the main steps 
involved in texture segmentation according to the 
invention, 

Fig. 3 is a flowchart to illustrate feature acquisition in 
a procedure according to Fig. 2, 

Fig. 4 is an illustration to determine a number of 
anisotropic scaling factors for an image point, 

Fig. 5 is a flowchart to illustrate texture classification 
in a procedure according to Fig. 2, 

Fig. 6 is an illustration to explain the initialization 
phase in a procedure according to Fig. 5, 
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Fig. 7 is an image sequence to illustrate the segmentation 
of Brodatz textures with an original image (a) , a 
filtered image (b) and a feature image (c) , 

Fig. 8 is an image sequence to illustrate the segmentation 
of a pixel image with two artificial textures with 
an original image (a), a depiction of the reference 
point for classification (b) , and a segmentation 
image (c) , 

Fi 9- 9 is another image sequence to illustrate the image 
g segmentation of an original image (a) with the use 

of varying distance measures (b, c) , 



in % 



Fig. 10 is an image sequence to illustrate texture 
segmentation in a grainy original image (a) and 
segmentations with varying distance measures (b-e) , 



fU Fig. 11 is a curve depiction to illustrate texture 

segmentation at varying levels of graininess, and 



Fig. 12 is a diagrammatic overview of a texture segmenting 
arrangement according to the invention. 

Texture segmentation according to the invention will be 
explained below using two-dimensional gray scale value 
images as an example, but is not limited thereto, but 
rather is applicable in the corresponding manner to any 
point distributions and feature combinations. For example, 
the point distribution can also be generated by several 
synchronously recorded time series of sensor signals, e.g., 
on a machine, wherein segmentation according to the 
invention is directed toward the search for specific time 
intervals within the time series, e.g., in which normal 
operating conditions of the machine or special faults are 
given. The examined point distributions can be continuous 
or discrete. The depicted image examples have been 



simplified in part for printing-related reasons, or 
provided with artificial structures (hatched lines or the 
line) , without constituting mandatory features of the 
invention. The individual steps involved in segmentation 
according to the invention will first be explained below 
drawing reference to Fig. 2 to 6. Subsequently, examples 
will be provided to illustrate the realization of the 
individual steps, and a device for implementing the 
procedure will be described. 



A) Segmentation Procedure 



A two-dimensional gray scale value image G(x,y) with size 
N • M is considered (N, M: number of image points or pixels 
in the x or y direction). A discrete gray value g(x, y) is 
allocated to each pixel (g e [0; 255]). The location and 
feature information allocated to each pixel generates a 
higher dimensional formation, which is a three-dimensional 
point distribution in this example. In this distribution, a 
three-dimensional vector p, = (x,y,g(x,y)) is allocated to 
each pixel. To ensure that the x-, y- and g-values lie in a 
comparable value range, it may be necessary to norm the 
gray values g ± . One possible norm takes the form of g norm = 
g • (N/255) . Hence, the pixel image is viewed as a point 
distribution in an artificial three-dimensional embedding 
space . 

Fig. la shows an example for a pixel image G. The pixel 
image with N = M = 256 contains four natural textures, 
which are readily identified upon visual inspection as a 
simple striped pattern (top right), regular honeycomb 
pattern (bottom left), completely irregular pattern with 
sharply delineated structures (top left) or irregular 
^indistinct" pattern (bottom right) . Image segmentation 
according to the invention is directed toward executing 
this texture recognition via the acquisition of local 
features of the pixel image for each image point, and 
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classifying the image points based on the acquired 
features. An overview of these steps is provided on Fig. 2 
with feature acquisition 100, texture classification 200 
and evaluation 300. Depending on the application, the 
evaluation may result in a further feature acquisition 
and/or texture classification. 

1) Feature Acquisition 

Feature acquisition 100 (see Fig. 3) is aimed at 
determining local features for each image point. The local 
features encompass characteristic image properties in the 
immediate environment of the image point, which is 
substantially smaller than the (global) overall image. 
Specifically, feature acquisition 100 involves determining 
the point distribution 110, determining scaling factors 
120, 130 and a feature extraction 140 to generate feature 
vectors, which are allocated to each image point. 

In the example in question, the determination of point 
distribution (step 110) involves a simple image recording 
and a known gray scale value evaluation (if necessary 
normed as specified) . In general, step 110 encompasses 
recording measured values based on the sensor technology 
selected as a function of application. 

With respect to determining the isotropic and anisotropic 
scaling factors (steps 120, 130), reference is made to the 
aforementioned DE-OS 43 17 746 and DE-PS 196 33 693. The 
procedures are known in the art, and will hence only be 
partially explained in detail here. 

First, the isotropic scaling factor (scaling index) is 
determined for each image point (step 120) . To this end, 
two spheres with differing radii ai, a 2 (ai < a 2 ) are 
concentrically placed around each point p, in the location 
and gray scale value space. A specific number of image 
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points are situated within each sphere, each with a gray 
scale value also designated as the overall mass M 
(respectively related to the sphere radius a lr 2 ) . 
According to equation (1), the isotropic scaling factor a 
is obtained as a logarithmic derivation of overall masses 
for both spheres: 

«(»i,yi|.„^)- loaW ''"- 1<> <'W-" (1, 

log a 2 - log a 1 



mit 

M(Xi, yi |a) = £©(a -||p. - p| 2 ), (2 ) 



wherein x±, y ± denote the coordinates of the image point in 
question, 0 the Heaviside function and |-| 2 the Euclidian 
norm. The scaling factor calculation represents a filter 
function, while the coordinates x ± , y ± denote the midpoint 
of the filter. 



Strong point density gradients around the image point in 
question yield a high scaling factor, while low gradients 
only produce a low scaling factor. The isotropic scaling 
factor a is here only characteristic for radial gradients, 
however. For applications with complex structures, it is 
necessary to also determine orientation properties for the 
local features. This is done by determining the anisotropic 
scaling factors (step 130) . 

The anisotropic scaling factors are determined similarly to 
the isotropic scaling factors, from gradients in the 
density of the environment points of a point in question, 
wherein the projections of the point numbers (projected 
masses M x , M y ) are ascertained from environment points to 
determine an orientation of the formation. To determine 
whether the environment points around an examined point 
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even have an orientation, it would essentially suffice to 
find two anisotropic scaling factors in a two-dimensional 
pixel image that each relate to the x or y axes of the 
pixel image- One characteristic of feature acquisition for 
image segmentation according to the invention involves 
ascertaining not just one value pair or value tuple from 
anisotropic scaling factors according to the dimensionality 
of the image in an examined image for each image point. 
According to the invention, several value pairs or value 
tuples of anisotropic scaling factors are determined 
according to the principles set forth with reference to 
Fig. 4. 

Fig, 4 shows an example of three differently oriented 
structures 41, 42 and 43 in a two-dimensional pixel image. 
The structures have a characteristic rotational angle <|) of 
0°, 45° or 90° relative to the x-axis. The determination of 
anisotropic scaling factors oe x or cc y is illustrated for the 
structure 41 by penciling in the spheres with radii a lr a 2 . 
The projection of environment points belonging to the 
structure 41 for the examined point p 4 i exhibit strong 
gradients in the y-direction, and weaker gradients in the 
x-direction. This results in a low scaling factor a x and a 
high scaling factor oc y . Conditions are correspondingly 
different for structures 42 and 43. 

It is shown that the anisotropic scaling factors can be 
used to show that the environment points have an 
orientation. However, the individual alignment of the 
orientation cannot be derived. This information can only be 
obtained if the anisotropic scaling factors for each 
examined point are determined in at least two coordinate 
systems rotated relative to each other. The scaling factors 
in the x, y and x*, y* coordinate systems (see structure 
42) yield additional information about the structural 
alignment given knowledge of the rotational angle (*90°) 
between the coordinate systems. 
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Therefore, the determination of anisotropic scaling factors 
in the procedure according to the invention (step 130) 
involves the acquisition of a scaling factor tuple 
comprised of several scaling factors each corresponding to 
the rotated coordinate systems for each examined image 
point. The principle of scaling factor determination in 
rotated coordinate systems shown in simplified form with 
reference to Fig. 4 is adapted to the application. This is 
represented on Fig. 3 by the step 131 of defining reference 
variables for scaling factor determination. These reference 
variables comprise the number of anisotropic scaling 
factors examined per image point, the sphere sizes and 
number and angle of coordinate system rotations. 

An angle scheme in which four rotational angles are set, 
each differing by 45°, has been proven advantageous for the 
evaluation of gray scale value images during coordinate 
system rotation. Therefore, one or two anisotropic scaling 
factors are determined for each image point in four 
coordinate systems each rotated by rotational angle <|>. It 
is emphasized that the number and level selected for the 
coordinate system rotations can be higher or lower 
depending on application. In this case, it does not matter 
that the respectively determined anisotropic scaling 
factors can be interpreted in a certain way or correlated 
with visually detectable image features. It is only 
important that several values be obtained for varying 
coordinate system rotations, since these values contain the 
complete information required for further texture 
classification (see below) . 

In order to calculate the anisotropic scaling factors, one 
proceeds from the projected local masses M x * according to 
equation (3) . 

M x *(Xi, yi |a, <|>) = ^ © (a - |x* - x*| 2 ) 0 (a 2 - ||p* -p*||). (3) 
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The axis represents the spatial direction within the 

rotated coordinate system to which the scaling factor is 
related. In equation (3), p* denotes the vectors of the 

image points in the rotated coordinate system. The 
transition from the original coordinate system (p - 
(x,y,g)) to the rotated coordinate system ( p* = (x*,y*,g*)) 
is achieved according to p* = D'p with the rotation matrix 
D from equation (4) . 



D = 



( cos <j> 
— sin <|) 



V 



0 



sin ( 

cos 
0 



0^ 
0 



(4) 



U 
n i 



Rotation matrix D is the known rotation matrix, wherein the 
g-axis is here the rotational axis. The second Heaviside 
function in equation (3) ensures that only points near the 
examined image point are taken into account. Similarly to 
equation (1), calculating the logarithmic derivation of 
projected masses M x * according to equation (5) yields the 
anisotropic scaling factors a: 



a <Xi, yi |,a lf a 2 ,4) = ^ M *« ^ »> I j^g ^ ¥ 

log a 2 - log a ± 



(5) 



The anisotropic scaling factor a according to equation (5) 
is calculated for each image point for four respectively 
different rotational angles 4. In the image examples 
described below, rotational angles 0°, 45°, 90° and 135° 
were used. Since the orientation information for each 
rotational angle is already contained in one of the two 
anisotropic scaling factors, which belong to a rotational 
angle, it is sufficient for subsequent texture 
classification if only one anisotropic scaling factor is 
determined for each rotational angle. 

In a preferred embodiment of the invention, the scaling 
factors were determined not just for a sphere pair, but for 
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two sphere pairs. In the 256 ' 256 pixel images examined 
here, a smaller scaling range with sphere radii ai = 2 and 
a 2 = 6 and a larger range with ai = 6 and a 2 = 12 proved to 
be advantageous. Additional scaling ranges (sphere pairs) 
can also be taken into account. 

Therefore, two isotropic and four anisotropic scaling 
factors that represent varying nonlinear filter functions 
for characterizing the textured image are obtained for an 
image point given two scaling ranges and four rotational 
angles. As demonstrated below, an outstanding segmentation 
result can already be derived using this set of ten scaling 
factors. Depending on the application, these parameters can 
be adapted to the specific task, or even optimized during 
the segmentation process, however. 

Feature acquisition 100 (see Fig. 3) is followed by feature 
extraction 140, which results in a feature vector x for 
each image point. The feature vectors x represent the 
input variables for the ensuing texture classification 200 
(see below) . 

In a simple case, the components of the feature vectors are 
formed by the local features determined for each image 
point. In the above example, the feature vector hence 
encompasses ten components comprised of two isotropic and 
eight anisotropic scaling factors. 

However, in a preferred embodiment of the invention, a 
statistical evaluation of the local features determined for 
the image points is first performed within the framework of 
step 141, followed by vector formation 142 out of local 
expected values. The accompanying local expectation value 
<ot> for each local feature is computed according to 
equation (6) . 
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< oc(xi,yi) 



> = ^Z a ( x 'y)©(|-|^ -x|)0(|-| yj -y|) (6) 



Parameter k here represents the size of a window shifted 
over the image points for taking neighboring local features 
into account. In the image examples described below, a 
window size of k = 40 was proven to be beneficial. As a 
result, ten expected values are formed for each image 
point, which take into account the local features of the 
adjacent image points. The advantage to this statistical 
processing of local features is that the ensuing texture 
classification is influenced significantly less by limiting 
effects on the image edges as a consequence. It has been 
shown that the limiting effects on the image edges need not 
even be considered at all to still achieve good 
segmentation results . 

Feature acquisition 100 ends with the determination of 
feature vectors x ± from the scaling factors or expected 
values for the scaling factors. Actual texture 
classification (clustering) 200 according to Fig. 5 now 
follows for additional image segmentation. 

2) Texture Classification (Clustering) 

The n feature vectors x ± form a set X according to 
equation (7 ) . 



In a first step 210, a finite number of texture classes 
arising in the image in question are first examined. These 
predetermined texture classes are respectively defined or 
entered from a storage medium, depending on application. An 
index i (i = 1, . . . , c) is allocated to each texture class 
(number c) . Depending on application, the texture classes 



X s {x x , x 



2' 



. . . , 




(7) 



are defined from experimental values, or within the 
framework of an optimization, in which image segmentation 
is performed repeatedly with varying class numbers and 
types . 

Within the framework of the partial supervision used 
according to the invention, it is now assumed that the 
allocation to a specific texture class is known for a 
specific, limited number of image points. Specific labels 1 
to c („labels M ) can hence be allocated to certain image 
points or image vectors. The remaining image points or 
image vectors are left without labels (^unlabeled") . 
According to equation (8), this can be written with 
reference to the set X introduced above: 



X nl' X l' 



,x 2 



nit 



r r 



' X nc' X 



1 ' 



(8) 



label 1 



label 2 



label c 



unlabeled 



X = 



U x u 



In equation (8), the superscript indices each relate to a 
label or the undesignated state without a label (u) . The 
subscript indices run from 1 to n (see equation (7)). The 
partial quantity X u is significantly greater than partial 
quantity X 1 of the feature vectors, for which the label is 
known . 



The allocation of reference points with labels to the known 
texture classes (step 220) is followed by the actual steps 
of the cluster procedure, namely the initialization phase 
230 and conversion phase 240. During the initialization 
phase 230, so-called ellipsoidal distance measures or a 
Euclidian distance measure is defined in the feature space, 
which establish specific scales for the reference points 
allocated to a respective texture class. During the ensuing 
conversion phase, the remaining image points are allocated 
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to the varying texture classes based on the distance 
measures or metrics. 

(2a) Initialization Phase 

The initialization phase 230 involves the steps of gravity 
center calculation 231, covariance matrix calculation 232, 
eigenvalue calculation 233, eigenvector calculation 234 and 
metrics definition 235. These steps are carried out in the 
high-dimensional feature space, which is fixed by the 
components of the feature vectors. In this example, the 
feature space is hence 10 dimensional. The steps involved 
in the initialization phase will be described below drawing 
reference to the simplified depiction in a two-dimensional 
feature space according to Fig. 6. 

In the feature space, the points of a texture each 
represent a texture class as cohesive formations, which are 
also referred to as clusters. Fig. 6 shows four clusters 
61-64 for two randomly selected components of the feature 
space according to the expected values for the scaling 
factors oci and <x 2 . The objective of the initialization 
phase 230 is to determine to which cluster, and hence to 
which texture class, an initially unclassified point 65 
(without label) is to be allocated. In a simple allocation 
procedure, a point without a label could simply be 
allocated to the cluster to which it has the lowest 
Euclidian distance. However, this can result in false 
allocations if the expansion and alignment of the cluster 
is not also taken into account. In the example shown, the 
point 65 has a slight distance to the cluster 62, but the 
cluster has a very characteristic longitudinal extension. 
For this reason, it might be more likely that point 65 
belongs to the more distant cluster 63, since this 
allocation is more compatible with the radial expansion of 
this cluster. In the initialization phase 230, a separate 
distance measure depending on the characteristic properties 



of the cluster orientation and formation is hence defined 
for each cluster, i.e., for each texture class. 

Gravity center calculation 231 takes place first. The 
cluster center is here calculated with the location vector 
u 1 in the feature space according to equation (9) 

-j ni 

o 1 = — £ x\ (9) 

for each cluster or each texture class i. The components 
(uj, u?;, ... u^) of the vector u 1 are the gravity center 

coordinates of the cluster formed by the feature vectors 

x k = ( x ik> x 2k>— x dk)- cl here denotes the number of 

dimensions in the feature space. The gravity center is 

symbolically entered on Fig. 6 for the cluster 61. 

Finally, step 232 involves calculating the covariance 
matrix C 1 for each cluster 1, whose elements C l rs are given 

by equation (10) . 

-1 ni 

C^ s = — Y (x x rk -o x r ).(x s k k -u*) with r, s = l, 2,..., d. (10) 

The elements of the covariance matrices establish a link 
between the deviation of components in each feature vector 
and the gravity center of the respective cluster i. In a 
simplified case according to Fig. 6, these constitute the 
deviations in the abscissa and ordinate directions. 
Matrices C 1 are symmetrical, thus enabling diagonalization 
and main axis transformation. 

In the ensuing step 233, the eigenvalues X\, A,* A,^ of the 

matrix C 1 are calculated for each cluster i. The square 
roots of the eigenvalues correspond to the standard 
deviations of the point distributions {x£}(k = 1, 2, 

n k ) relative to the respective main axes. Finally, the 
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eigenvectors for each matrix C 1 are calculated in step 234. 
The eigenvectors form the matrices D 1 according to equation 
(11) . 



D i = 
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Matrices D 1 describe the transition from the original 
coordinate system shared by all clusters to a cluster- 
specific coordinate system fixed by the main axes of the 
respective cluster i. This coordinate system is illustrated 
based on the cluster 64 on Fig. 6, for example. 

Therefore, a local coordinate system whose axes are 
calibrated by the formation of the respective cluster is 
introduced for each texture class introduced within the 
framework of partial supervision in the initialization 
phase 230. The local calibrations yield information about 
the orientation and formation of the respective clusters, 
and hence make it possible to define cluster-specific 
distance measures local to the feature space (ellipsoidal 
metrics) . These distance measures are defined during step 
235. The coordinates of the cluster centers, the square 
roots of the eigenvalues and the eigenvectors are stored 
for each texture class or each cluster. This yields N var = i 
* (2 • d + d 2 ) variables according to the number of various 
clusters (classes) i and the dimension of the feature space 
d, for each cluster d, since d parameters X±, d location 
coordinates of the cluster centers and d 2 eigenvectors are 
taken into account for each cluster. These variables 
determine the ellipsoidal distance measures (metrics) in 
the feature space, in which the cluster allocation of the 
remaining points without labels takes place (see also 
equation (15) ) . 
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It must be noted that the distance measure definition with 
steps 232 to 234 covered herein is not a compulsory feature 
of the invention. Image segmentation based on the feature 
vectors determined via nonlinear filtering during feature 
acquisition 100 can also be realized using a simpler 
distance measure, e.g., based on the Euclidian distance 
between an image point and the gravity center of a cluster. 
Accordingly, steps 232 to 234 could be skipped (dashed 
arrow on Fig. 5) . However, the selection of distance 
measure impacts the quality of image segmentation (see Fig. 
10, 11) . 

(2b) Conversion Phase 

In the conversion phase 240, previously unclassified image 
points are allocated to a cluster based on the distance 
measure belonging to this cluster. To this end, the 
distance vectors y^* for each of the y cluster centers are 

calculated in step 241 for each image point without a label 
x^ (1 = 1, 2, n u ) according to equation (12) . 

yf = - o 1 . (12) 

The distance vectors y\* are subsequently transformed into 
the coordinate systems of the main axes for each cluster 
(step 242) . This is done according to equation (13) using 
the transfer matrices D 1 . 

^ = D 1 • yf . (13) 

As a result, c distance vectors corresponding to the c 
local coordinate systems of the clusters are present for 
each unclassified image point. In each coordinate system, 
the components of the vectors yi* are calibrated relative 
to the standard deviations of each cluster using the 



eigenvalues of the respective local covariance matrix C 1 
(step 243) . 



This calibration is performed according to equation (14) : 



Based on the calibrated components, the distances of each 
image point to be classified to each of the i cluster 
centers are calculated according to equation (15) (step 
244) . 



Finally, the feature vector x* (or accompanying image 

point) is allocated to the cluster or texture class for 
which variable A\ has the lowest value (step 245) . As a 

result, the number of all image points is completely 
classified, i.e., each image point is allocated to a 
cluster or a texture class. 

If classification takes place with a simple Euclidian 
distance, step 241 encompasses calculating the vectors of 
an image point for each gravity center of the clusters. 
Then, one can proceed directly to step 241. 

(3) Evaluation 

Within the framework of subsequent evaluation 300 (see Fig. 
2) , there now takes place a segmentation step in which 
image points that have been allocated to a shared cluster 
are provided with a shared label, stored together and/or 
marked in a display of the processed image by a false color 
depiction. Simultaneous to this segmentation, the image 
points belonging to a texture class also readily permit a 




(14) 




(15) 
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size determination of the partial area of the image formed 
by a texture according to one of the known numerical 
evaluation procedures . Additional steps during evaluation 
300 include the relaying of texture data to various 
auxiliary processors and/or displays configured as a 
function of application. A partial area can be self- 
contained or consist of several segments separated from 
each other. 

Compliance with a quality standard can also be checked 
within the framework of the evaluation 300 , and given a 
negative result, a return to feature acquisition 100 or 
texture classification 200 can be effected with the 
establishment of new parameters. Quality standards include 
correctly classified image points known in advance. To this 
end, the procedure outlined above can be modified so that 
not all known image points are taken into account during 
step 220 and incorporated in to the ensuing initialization 
phase 230. Initialization can be performed with a first 
portion of known image points with a known texture 
classification. The result of the conversion phase 240 can 
then be checked using the remaining known image points. 

B) Experimental Results 

Results of image segmentation according to the invention 
based on the example of gray scale value images with 
natural or artificial textures will be presented below 
drawing reference to Fig. 1 and 7 to 11. Fig. 1 shows a 
256 * 256 pixel image with four natural textures that were 
already specified above (Fig. la) . These involve the 
Brodatz textures D03, D90, D93 and D51. Fig. lb shows the 
segmentation result. 97.3 % of the image points were 
correctly classified. As also revealed, the boundaries 
between the textures are reproduced relatively well. 
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Fig. 7 shows additional Brodatz textures in a 128 * 256 
pixel image. The original textures D55, D68 are depicted on 
Fig, 7a, Fig. 7b shows the image of the anisotropic scaling 
factors filtered with the parameters ai =6, a 2 = 12 and § 
= 45°. Fig. 7c is obtained after a smoothing according to 
the aforementioned step 141. 

The image sequence on Fig. 8 illustrates the image 
segmentation according to the invention based on the 
example of artificial textures. Fig. 8a shows two 
artificial textures, of which the first texture in the left 
half of the image consists of triangles, and the second 
texture in the right corner of the image consists of 
arrows. The textures have the same second-order statistics, 
and can therefore not be distinguished between using local 
linear feature acquisitions (see B. Julesz in „Rev. Mod. 
Phys. x \ Vol. 63, 1991, p. 735). The white image points 
highlighted for the sake of clarity on Fig. 8b illustrate 
the reference points or image points with labels used for 
initializing the cluster procedure. These are 655 pixels 
selected completely at random. The share of reference 
points hence lies at 1 % relative to the overall point 
number (these conditions were also used for the remaining 
examples) . Finally, Fig. 8c shows that, despite the 
identical second-order statistics, reliable image 
segmentation can be achieved using the procedure according 
to the invention. 98.2 % of the image points have been 
correctly classified. 

Fig. 9a shows another set of four natural Brodatz textures 
(D96, D55, D38 and D77) . Fig. 9b and 9c illustrate the 
different results when using various distance measures for 
texture classification. When using ellipsoidal metrics, 
97.1 % of the points are correctly classified according to 
Fig. 9b. By contrast, if only a Euclidian metric is used 
for texture classification, only 94.9 % of the image points 
are correctly classified according to Fig. 9c. As a result, 



the segmentation result is worse in comparison to Fig. 9b, 
but still adequate for various applications. 

If the texture image is now overlaid with noise, Fig. 10a 
results. The noise is an additive white Gaussian noise with 
a noise level corresponding to a predetermined standard 
deviation a nC ise- In the depiction according to Fig. 10a, a 
noise = 0.1 • AI, wherein AI denotes the intensity range of 
the undisturbed image. 

When using ellipsoidal metrics according to Fig. 10b, an 
acceptable segmentation result is still achieved for the 
grainy image, while the result is heavily skewed when using 
a Euclidian metric according to Fig. 10c. Given an increase 
in the noise percentage to a no ise = 0.2 * AI, there is a high 
stability in image segmentation, provided it is based on 
the ellipsoidal metrics (Fig. lOd) . Image segmentation 
based on the Euclidian metric exhibits a further 
deterioration, however (Fig. lOe) . 

The noise effect visualized on Fig. 10 is also illustrated 
in the curve depiction according to Fig. 11. Fig. 11 shows 
the number of correctly classified image points (ccp) in 
percent as a function of the noise level a n0 ise- The solid 
line corresponds to the segmentation result when using the 
ellipsoidal metrics, while the dashed line reflects the 
result obtained when using the Euclidian metric. The 
significantly higher stability of the image segmentation in 
the first case is visible. 

The image examples highlight the following significant 
advantages to image segmentation according to the 
invention. First, it must be pointed out that all texture 
recognitions are based only on knowledge of the feature 
image (feature vectors) . In image segmentation, the 
specific image point coordinates are not the critical 
factor, but only the properties of the feature vectors. The 



procedure is parameter-free relative to cluster allocation. 
This represents a significant difference in comparison to 
the parameters to be optimized in conventional texture 
classification procedures. 

C) Segmentation Arrangement 

The diagrammatically illustrated segmentation arrangement 
on Fig. 12 consists specifically of a measuring device 1, 
with which the point distribution to be segmented is 
acquired, a filtering device 2, which is designed to 
evaluate the measured point distribution according to the 
filtering functions explained above, and to this end 
comprises scanners, counters and means for acquiring or 
statistically processing scaling factors, an input device 
3, which is used to enter information for the partial 
supervision of image segmentation, i.e., predetermined 
textures and the group of reference points with labels, a 
calculating device 4 for determining and evaluating 
distance measures for the texture features of all remaining 
points, and an output device 5 to display, temporarily 
store or relay texture segments. Additional components of a 
device according to the invention, e.g., a controller, are 
not shown. The filtering device 2 and calculating device 4 
are preferably formed by a shared computer. 
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CLAIMS 

1. A method for segmenting a point distribution with 
numerous points into partial areas, which each exhibit 
specific structural elements, with the following 
steps : 

a) for each point (Pi), a feature vector (x^) is 
determined, whose components are determined based 
on several scaling factors belonging to the 
respective point; 

b) the accompanying feature vectors (p^;) are determined 

for a predetermined number of reference points (xj) 
of the point distribution for which the allocation 
to one of the structural elements is given, and 
texture classes each corresponding to the 
underlying structural elements are formed out of 
the feature vectors of the reference points; 

c) for all remaining points (Pi) of the point 
distribution that are not reference points, the 
distance between the respective point and each of 
the texture classes is determined based on a 
distance measure in the feature space, which is 
fixed by the components of the feature vectors; 

d) the points (p£) are each allocated to the texture 
class for which the lowest distance was determined; 
and 

e) the partial areas of segmentation are formed out of 
the respective reference points belonging to a 
texture class and the points allocated in step d) . 



2. A method according to claim 1, in which several 
isotropic and anisotropic scaling factors (a) are 
determined as components of the feature vectors (5^). 



A method according to claim 1, in which expected 
values of several scaling factors are determined as 
components of the feature vectors (x^. 

A method according to claim 2 or 3, in which the 
anisotropic scaling factors ascertained for a point 
are determined in various coordinate systems rotated 
relative to each other. 

A method according to one of the preceding claims, in 
which an ellipsoidal distance measure is used as a 
local distance measure in step c) for each texture 
class . 

A method according to one of claims 1 to 5, in which a 
Euclidian distance measure is used as a shared 
distance measure in step c) for all texture classes. 

A method according to one of the preceding claims , in 
which the segmented partial areas are displayed, 
temporarily stored and/or further processed. 

A method according to claim 7, in which the size of 
the partial areas is quantitatively acquired 
simultaneous to displaying the segmented partial 
areas . 

A image segmentation arrangement for segmenting a 
point distribution out numerous points, with a 
measuring device (1) for acquiring the point 
distribution, a filtering device (2) for scanning and 
filtering the point distribution, an input device (3) , 
a calculating device (4) and an output device (5), 
wherein the filtering, input and calculating devices 
(2, 3, 4) are designed to execute a method according 
to one of claims 1 to 8 . 



A use of a method or a device according to one of the 
preceding claims to process: 

- images of medical or biological objects; 

- images of materials; 

- point distributions of complex static systems; 

- point distributions that reproduce system status of 
complex, dynamic systems; and 

- time patterns for dynamic systems. 
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